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Decays of the ^/'(2S) into Vector plus Pseudoscalar meson final states have been studied with 14 
million '4>{2S) events collected with the BESII detector. Branching fractions of ^{2S) (j>Vi 4'V\ s-nd 
uri', and upper limits of (25) cpn'^ and wr? are obtained: B{'^{2S) (prj) = (3.3±1.1±0.5) x 10"^, 
B(ij{2S) 4>'n') = (3.1 ± 1.4 ± 0.7) X 10-^ and B{iP{2S) ojt]') = (3.2t|^ ± 0.7) x 10"^ and 
B{ip{2S) -y 07r°) < 0.40 x 10"'"', and B{ip{2S) ^ aJT]) < 3.1 x 10"^ at the 90 % C.L.. These results 
are used to test the pQCD "12% rule". 

PACS numbers: 13.25.Gv, 12.38.Qk,14.40.Gx 



I. INTRODUCTION 

It is expected in perturbative QCD that both J/tp 
and V'(2'S') decays to light hadrons proceed dominantly 
via three gluons or a single virtual photon, with widths 
proportional to the squares of the cc wave functions at 
the origin [1], which are well determined from dilepton 
decays [2]. This led to the "12% rule", i.e. 

_ E(V>(2g) ^ h) ^ B{i,{2S) ^ e+e-) ^ 
^''^ B{J/^^h) - B{J/^^e+e-) 

A strong violation of this conjecture was first observed 
by the Markll experiment in the Vcctor-Pseudoscalar 
meson (VP) final states, pn and K*+{S92)K- [3]. 
Significant suppressions observed in four Vector- 
Tensor decay modes (a;/2(1270), pa2(1320), (pf2, and 
K*{892)K*{1430) + cc.) [4] make the puzzle even 
more mysterious. Numerous theoretical explanations 
have been suggested [5], but the puzzle still remains 
one of the most intriguing questions in charmonium 
physics. Recently both CLEO and BES reported new 
measurements of VP channels [6-8] with higher statis- 
tics and confirmed the severe suppression for pn and 
K*{892)K + C.C.. 

In this letter, we report measurements of tlj{2S) de- 
cays into 5 VP channels: (^7r°, (prj, (prj' , urj, and urj' 
using 14 million 'ip{2S) events collected with the BE- 
SII detector, where the branching fractions of (prf and 
cur]' are the first observations. The results are com- 
pared with the corresponding J/ ^ branching fractions 
to test the "12% rule". Also, the branching fractions 
provide useful informations on DOZI suppressed de- 
cay ip{2S) 07r° and on the quark components of rj 
and r]' [9]. 

II. THE BESII DETECTOR 

The Beijing Spectrometer (BESII) is a conventional 
cylindrical magnetic detector that is described in de- 
tail in Ref. [10]. A 12-layer Vertex Chamber (VC) 
surrounding the beryllium beam pipe provides input 
to the event trigger, as well as coordinate informa- 
tion. A forty-layer main drift chamber (MDC) located 



just outside the VC yields precise measurements of 
charged particle trajectories with a solid angle cov- 
erage of 85% of 47r; it also provides ionization energy 
loss (dE/dx) measurements which are used for particle 
identification. Momentum resolution of 1.7%a/1 + p"^ 
{p in GeV/c) and dE/dx resolution for hadron tracks 
of 8% are obtained. An array of 48 scintillation 
counters surrounding the MDC measures the time of 
flight (TOP) of charged particles with a resolution of 
about 200 ps for hadrons. Outside the TOF coun- 
ters, a 12 radiation length, lead-gas barrel shower 
counter (BSC), operating in limited streamer mode, 
measures the energies of electrons and photons over 
80% of the total solid angle with an energy resolution 
of ge/E = 0.22/v^ {E in GeV). A solenoidal mag- 
net outside the BSC provides a 0.4 T magnetic field 
in the central tracking region of the detector. Three 
double-layer muon counters instrument the magnet 
flux return and serve to identify muons with momen- 
tum greater than 500 MeV/c. They cover 68% of the 
total solid angle. 

In this analysis, a GEANT3 based Monte Carlo 
package with detailed consideration of the detector 
performance (such as dead electronic channels) is 
used. The consistency between data and Monte Carlo 
has been carefully checked in many high purity physics 
channels, and the agreement is reasonable [11]. 

III. EVENT SELECTION 

The data sample used for this analysis consists of 
(14.0 ± 0.6) x lO'^ '4){2S) events [12], collected with 
BESII at the center-of-mass energy ^/s = m,^(2S) ■ The 
decay channels investigated are ip{2S) into 07r°, (j)r], 
(j)r]' , u)i], and our]', where (p decays to K^K~, lo to 
7r'''7r~7r°, rj' to rjn'^TT' or 77r+7r~, and tt" and rj to 2j. 
The events have either two or four charged tracks plus 
n (n > 1) photons. 

A neutral cluster is considered to be a photon can- 
didate if the following requirements are satisfied: it 
is located within the BSC fiducial region, the energy 
deposited in the BSC is greater than 50 MeV, the first 
hit appears in the first 6 radiation lengths, the angle 
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I cos^l < 0.8. Here Vx, Vy, and Vz are the x, y, and 
z coordinates of the point of closest approach of the 
track to the beam axis. 

The TOF and dE/dx measurements for each 
charged track are used to calculate Xpioi'^) values 
and the corresponding confidence levels Probpioii) 
for the hypotheses that a track is a pioii. kaon, or 
proton, where i {i = 7:/K/p) is the particle type. 
For events with K^K~ decays, charged kaon 
candidates are required to have ProbpioiK) larger 
than 0.01 or larger than Probpujlir) and Prohpiij{p)] 
while for events with lo — » 7r+7r~7r° decays, at least 
half of the charged pion candidates in each event are 
required to have ProbpioiTr) > 0.01. 



FIG. 1: Scatter plot of rriKK versus m^-y for ^p{2S) 
is"*" if ~ 77 candidate events. The two ellipses in the plot 
indicate the 2a contours of the signal regions. 
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FIG. 2: Distribution of 77 invariant mass for '(p{2S) 
077 candidate events (solid histogram) and Monte Carlo 
simulation of ip{2S) — » (jm" (dotted histogram) with ar- 
bitrary normalization. The curve shows the best fit de- 
scribed in the text. 

in the xy plane (perpendicular to the beam direction) 
between the cluster and the nearest charged track is 
greater than 16° (this requirement is not applied for 
channels involving more than two photons), and the 
angle between the cluster development direction in the 
BSC and the photon emission direction from the beam 
interaction point (IP) is less than 37°. 

Each charged track is required to be well fitted by a 
three-dimensional helix, to have a momentum trans- 
verse to the beam direction greater than 70 McV/c, 
to originate from the IP region, Vxy = \J^x + < 2 
cm and \Vz\ < 20 cm, and to have a polar angle 



A. ip{2S) <j)n° and </)T7 

The K^K~'yj final state is utilized to measure 
these two channels. Two good charged tracks with 
net charge zero and at least two photon candidates 
are required. Next a four constraint (4C) kinematic 
fit (xline) to the K^K~'yj hypothesis is performed, 
and the confidence level of the fit is required to be 
larger than 0.01. If there are more than two pho- 
tons, the fit is repeated using all permutations of pho- 
tons, and the two photon combination with the min- 
imum Xkine selected. This procedure is used for 
all channels. To suppress backgrounds from ir/K/p 
misidentification, the combined chisquarc [13], xtomb^ 
for the tp{2S) K^K~jj assignment is required 
to be smaller than those of iIj{2S) — > 7r+7r~77 and 
tp{2S) — > ppjj- Figure 1 shows the scatter plot of the 
invariant mass of K~^K~ (fTT-K+K-) versus that of the 
two gammas {m-yy). A clear cluster can be observed 
in the (pr] signal region, while only one event appears 
in the (pw^ region. 

To select cfi decay candidates, wc require \m,j^+j^- — 
1.02| < 0.02 GeV/c^. Here the experiment mass res- 
olution for (j) is 2.5 MeV/c^. The 77 invariant mass 
distribution for events with (p candidates is shown in 
Fig. 2. Fitting this distribution with tt'^ and rj bmc- 
tions determined from Monte Carlo simulation, plus 
an r]K^K~ background determined by the sideband 
and a first order polynomial to describe phase space 
backgroimd, 16.7 ± 5.6 (fyrj events arc obtained with 
the statistical significance of 3.6(7 [14]. While for the 
(jm^ channel, the observed events and the estimated 
backgromid in the signal region are 4 and 6.2, respec- 
tively, which corresponds to the upper limit of 4.4 ^tt" 
events at the 90% confidence level [15]. 
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FIG. 3: The 77r+7r~ invariant mass distribution for 
ip{2S) — > 07r'''7r~7 candidate events. The curve shows the 
best fit described in the text. 
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FIG. 4: The -k^-k't] invariant mass distribution for 
ip{2S) 4>n~^n~ri candidate events. The curve shows the 
best fit described in the text. 

B. ■0(25) <t>n' 

Two decay modes of r?'(958) are used, r?'(958) 
tt+tt"?/ and 777 "'"7r~. Their final states are 
K^K~-K^'K~ + m7, where m = 1 for r/' — > 77r+7r~ 
and TO = 2 for 77' — > tt+tt"/;, 77 — > 77. Events 
with four charged tracks with net charge zero and at 
least m photon candidates are selected. A 4C kine- 
matic fit is performed for the hypothesis tp{2S) 

K^TT'^TT~mj, and the confidence level of the fit 
is required to be larger than 0.01. The correspond- 
ing Xcomb fo'^ the V(2'S') — > K~^K~ Tr~^Tr~m'y assign- 
ment is required to be smaller than those of tp{2S) — > 
7r^7r^7r+7r^m7 and iJj{2S) K^m"f. 

The additional requirement \mx+K- — 1-02| < 0.02 
GeV/c^, is used to select ^ candidates, and the 
TO^^+^- spectrum of selected events is shown in Fig. 3. 
By fitting this spectrum with an rj' function, plus 
an rj'K^K^ background determined by the sideband 




FIG. 5: The n^Tr~n invariant mass distribution for 
7r''"7r~7r°r? candidate events (solid histogram) and Monte 
Carlo simulation for 'tp{2S) tor] (dotted histogram) with 
arbitrary normalization. The curve shows the best fit de- 
scribed in the text. 

and a second order polynomial for phase space back- 
ground, 5.8 ± 3.2 (f)ri' candidate events are obtained. 
The statistical significance is about 2.0(t. Here, the 
shape of rj' is determined from Monte Carlo simula- 
tion of il>{2S) — > ^T]', (j) — > K^K~ , and 77' — > 77r+7r~. 

For the final state isr+i^~7r+7r~77, the 
K~ -K^ 'iT~ r] candidate events are required to 
satisfy |to^^ - 0.547| < 0.05 GeV/c^. Back- 
grounds from tjj{2S) — > 7r+7r~J/V', J/ip — > ^r? and 
■i/'(25) vJ/'^'j J/i' ^ (/)7r+7r~ are eliminated with 
two additional requirements |m^ecoi( ~ 3.1| > 0.1 
GeV/c^ and mj^+j^-^+^- < 2.9 GeV/c^, respectively. 

With the requirement \mK+K- — 1-02| < 0.02 
GeV/c^, the 7r+7r~77 invariant mass distribution for 
<pTr^Tr~r] candidate events is shown in Figure 4. Fit- 
ting the spectrum with an r]' function plus a second 
order polynomial for background, 2.6 ± 1.8 ^77' candi- 
date events are obtained. The statistical significance 
is about 2.0(7. The rj' shape is determined from Monte 
Carlo simulation of iIj{2S) cjn]' , (j) K^K~ and 
T}' — > 'iT~^iT~r], 7] — > 77. 

C. ip{2S) WT) 

Here, the final state studied is 7r+7r~7777. Events 
with two charged tracks with net charge zero and four 
or five photon candidates are selected. A 4C kine- 
matic fit is performed for the hypothesis iJj{2S) — > 
7]-+7r^7777, and the fit confidence level is required 
to be larger than 0.01. To remove background from 
tt/K misidentification, xlomb required to be smaller 
than that of iJj{2S) K~jjjj. Candidate events 

must satisfy Im^j-yj — 0.135| < 0.05 GeV/c^ and 
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FIG. 6: The 7p° invariant mass distribution for '4>{2S) — > 
Tr~^ ■K~ tt'^ ■K~ ^ candidate events. The curve shows the 
best fit described in the text. 

I'^7374 ^ 0.547| < 0.05 GeV/c^ for the four photon 
candidates (71727374), where the subscripts permute 
over all six combinations. Events with one and only 
one combination satisfying the above criteria are kept 
for further analysis. 

Figure 5 shows the 7r"'"7r~7r°(TO,rOT7r) invariant mass 
distribution after the above selection; no clear uir] sig- 
nal is seen. The distribution is fitted with an oj signal 
determined by Monte Carlo simulation and a poly- 
nomial background, the observed events and the esti- 
mated background in the signal region are 23 and 25.0, 
respectively, which corresponds to the upper limit of 
9.7 wjj events at the 90% confidence level. 

D. %I>{2S) wt;' 

Two 77' (958) decay modes are used in this mea- 
surement, similar to the measurement of ^^{23) — > 
(j)?]' . Final states studied are 7r+7r~7r"''7r~ + rwy, 
where m = 3 for rj' ^'!t~^tt~ and to = 4 for 
7]' — > ?77r+7r~. Events with four charged tracks with 
net charge zero and m. or m, + 1 photons candidates 
are selected. A 4C kinematic fit to the hypothesis 
ip{2S) 7r+7r~7r+7r~TO7 is performed, and its con- 
fidence level is required to be larger than 0.01 and 
larger than that of tp{2S) K^K~TT'^TT~m/y to sup- 
press possible backgrounds due to particle misidenti- 
fication. Backgrounds from ip{2S) — > tt+tt" J/V' are 
rejected with the requirement that the mass recoiling 
from every tt+tt" pair satisfy IWrecoi; ~ "^•-'-1 ^ ^-^^ 
GeV/c^. 

For 7r+7r~7r+7r~777, one and only one pair among 
the three good photon candidates is required to satisfy 
|to-y^ — 0.135| < 0.05 GcV/c^: this pair is taken as a 
7r°. To avoid contamination from ip{2S) cott^tt^, 
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FIG. 7: The r]-iT~^-iT~ invariant mass distribution for 
tl>{2S) — > n'^Tr~n'^Tr~n°ri candidate events (solid his- 
togram) and Monte Carlo simulation of ip{2S) — > 
ujr) ,r] —^ 7v'^TT~ri (dotted histogram) with arbitrary nor- 
malization. The curve shows the best fit described in the 
text. 



the 7r+7r^7r'^7r+7r^ invariant mass should be less than 
3.5 GeV. Since the dominant decay of 77 into 771 "'"tt" is 
7/9, an additional requirement |m^+^- — 0.771| < 0.15 
GeV is applied to select 7r+7r~7r''7p candidates, where 
ttJ^tt^ is any combination from the four charged pion 
candidates. 

An additional requirement \m„+^-T^o — 0.7826| < 
0.05 GeV /(? is used to select w candidates. The 7p° 
invariant mass spectrum for selected events is shown 
in Fig. 6. It is fitted with an 77' function determined 
by Monte Carlo simulation for tp{2S) uij]', 1]' 
7/9, plus a second order polynomial for background, 
as shown in Fig. 6; 4.2 ± 2.7 events are obtained with 
a statistical significance of 1.9cr. 

For the final state 7r"'"7r~7r"'"7r~7777, the selection 
of 7r° and r] is the same as for the unj channel. Only 
events with only one tt candidate (| 0.1351 < 

0.05 GeV/c^) and one 77 candidate {\mj^ - 0.547| < 
0.05 GeV/c^) from amongst the four photons are kept 
for further analysis. 

An additional requirement \'ni^+^- ^0 — 0.78261 < 
0.05 GeV is made to select lu candidates. The 7y7r+7r~ 
invariant mass spectrum is shown in Fig. 7. It is fitted 
with an rj' function determined by Monte Carlo sim- 
ulation for ^j(2S) Lurj' , rj' rjTT^Tr' , plus a second 
order polynomial for background, as shown in Figure 
7. With 1 event observed in the signal region and 
3.2 background events estimated from sidebands, the 
candidate curj' signal is Ota'' event (at 68.3% C.L.) 
assuming the Poisson variate [15]. 
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IV. SYSTEMATIC ERRORS 

Many sources of systematic error are considered. 
Systematic errors associated with the efficiency are 
determined by comparing J/ip and ip(2S) data and 
Monte Carlo simulation for very clean decay channels, 
such as tp{2S) tt+tt" J/V', which allows the determi- 
nation of systematic errors associated with the MDC 
tracking, kinematic fitting, particle identification, and 
photon selection efficiencies [4, 11]. 

Since the decay rj' to 77r+7r~ includes -yp^ and 
non-resonant 77r+7r~, the tt+tt" invariant mass spec- 
trum in the Monte Carlo simulation is obtained from 
J/tjj ^ 77?', T]' — > 77r+7r~ data. The uncertainty of 
their detection efficiency from tt~^tt~ invariant mass 
spectrum is 3%, which is included in systematic er- 
rors. 

The uncertainties of the branching fractions of in- 
termediate states, the background shapes in fitting, 
and the total number of Tp{2S) events are also sources 
of systematic errors. Table I summarizes the system- 
atic errors for all channels. 

Contributions from the continuum e~^e~ ^ 7* — > 
hadrons [16] are estimated using a data sample of 
6.42 ± 0.24 pb-i taken at ^/s = 3.65 GeV [17], cor- 
responding to about one-third of the integrated lumi- 
nosity at the ^p{2S). Since no signal is observed for 
any channel analyzed, the continuum contribution and 
possible interference are not taken into consideration. 

TABLE I: Summary of relative systematic errors (%). 
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1.9 


Background shape 


0.0 


11.0 


16. 


0.0 


18.9 


MC statistics 


1.4 


2.1 


1.8 


1.3 


2.0 


Branching ratios 


1.4 


1.6 


3.8 


1.0 


3.5 




4.0 


Total 


9.6 


15. 


21. 


11. 


23. 



where X is the intermediate state, Y the final state, 
and e^'-' the detection efficiency. 

Table II summarizes the observed numbers of 
events, detection efficiencies, and branching fractions 
or upper limits for the channels studied. The branch- 
ing fractions of il>{2S) — > ^r]' and ^^{28) — > wq' are 
calculated from the sum of events observed in the 
rj' 7/3 and r/TT+Tr" channels , and an efficiency de- 
termined from the individual efficiencies weighted by 
the branching fractions of these two channels. The 
upper limit for ip{2S) Lurj' branching fraction at 
90% confidence level is 9.2 x 10~^. For comparison, 
the table includes the corresponding branching frac- 
tions of J/tp decays [2], as well as the ratios of the 
ip{2S) to J/tp branching fractions. Decays of ip{2S) 
to (j)r] and cor] are suppressed by a factor of 2 and 12, 
respectively, compared with the 12% rule, while (prj' 
and LOT]' are consistent with the rule within large er- 
rors. It is worth pointing out that the ratio of gj^^^)) 
is 2.3 ± 0.3 and 1.1 ± 0.7 for J/ip and tp{2S) decays, 
respectively, which are consistent within 2a, while the 
ratio of ^^^^^}^ is 9. 5 ±1.7 for J/i/; decay, which is much 
larger than that of tp{2S) decay. 

This analysis without considering the continuum 
contribution (although not seen at present measure- 
ment) and possible interference might bring some un- 
certainty, which could only be clarified later by more 
accurate experiments such as CLEO-c or BESIII [18]. 
If the continuum contribution is treated incoherently, 
the continuum events, assuming the Poisson distribu- 
tion, for (pri, (prj' and wry' channels are Ot^-° at 68.3% 
confidence level with the normalized integrated lumi- 
nosity, this yields the branching fractions of ■ip{2S) 
(pT], (prj' and ujtj' to be (3.3;i;4±0.5) x 10"^, (3.ll2:o± 
0.7) X 10-5 and (3.2^2-4 ± o.7) x 10"^, respectively. 

In conclusion, the branching fractions for tp{2S) — > 
(prj, cpT]' , and oji]' and upper limits for (pn'^ and ojr] 
are presented. Our results for (pij' and cor]' are first 
measurements, while our results for ^tt^, (pr], and oor] 
are consistent with those of CLEO [6]. 
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TABLE 11: Branching fractions and upper limits (90% C.L.) measured for ip{2S) — > Vector + Pseudoscalar. Results for 
corresponding J/ip branching fractions [2] and the ratios Qh ~ ^sn^f^j^h) given. 



h 




e 


B{xP{2S) -*) 
xlO"^ 


B{J/%1} -*) 


Qh 






(%) 


xlO-" 


(%) 


dm" 


< 4.4 


16.1 


< 0.40 


< 0.068 






16.7 ±5.6 


18.9 


3.3 ± 1.1 ±0.5 


6.5 ±0.7 


5.1 ± 1.9 


(f>ri'{r]' — > 777 "'"Tr") 


5.8 ±3.2 


11.1 








4>r)'(ri' —^ r)'K'^n~) 


2.6 ± 1.8 


3.8 








(fn}' {combined) 


8.4 ±3.7 


8.4 


3.1 ± 1.4 ±0.7 


3.3 ±0.4 


9.4 ±4.8 


ujr) 


< 9.7 


6.3 


< 3.1 


15.8 ± 1.6 


< 2.0 


ijjvl (rj — > 77r"'"7r~) 


4.2 ± 2.7 


2.6 








0177' (r?' — > r77r"'"7r~) 


n n+^-^ 


1.8 








ujr]' {combined) 


4 2+^-2 

^.Z_2 7 


2.3 


3.2l2;o ±0.7 


1.67 ±0.25 
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